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The mass and wave function of a six-quark system with quantum numbers 
J p =0", T=0, called d', are calculated. We use a colored diquark-tetraquark 
cluster model for the six-quark wave function. A constituent quark model 
Hamiltonian with a two-body confinement potential, and residual one-gluon, 
one-pion, and one-sigma exchange interactions is used. The complications due 
to the quark exchange interactions between tetraquark and diquark clusters 
(Pauli principle) are taken into account within the framework of the Res- 
onating Group Method. The calculated d! mass is some 350 MeV above the 
empirical value if the same two-body confinement strength as in the nucleon 
and A is used. This paper also examines the validity of the usual assumption 
of a universal two-quark confinement strength. We propose that the effective 
two-body confinement strength in an exotic six-quark system, such as the d' , 
could be weaker than in a single baryon. The weaker confinement hypothesis 
leads to a d' mass of = 2092 MeV and a d! radius of r> = 1.53 fm. 
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I. INTRODUCTION 



First indications for the existence of the d! dibaryon came from the narrow peak observed 
in the pionic (7r + ,7r~) double charge exchange (DCX) reaction on nuclei. At an incident pion 
energy of T n ~ 50 MeV and forward pion scattering angles = 5°, the DCX cross section, 
for example for tt + + 12 C — > 12 + tt~ , displays a narrow peak |Ih|3|. Because of 
charge conservation, the DCX reaction involves at least two nucleons within the nucleus, 
and it has been shown that the DCX reaction is therefore very sensitive to short-range NN- 
correlations 0]. The possible importance of explicit quark degrees of freedom in DCX at 50 
MeV and forward angles was pointed out some time ago ||. In the meantime, dedicated 
DCX experiments on a large number of light and medium heavy nuclei ranging from 

7 Li to 56 Fe have unambiguously confirmed the existence of a narrow resonance-like structure 
at this pion energy and angle. 

While conventional DCX calculations |5]|| have so far been unable to explain these 
experimental results, the assumption of a single narrow baryon number B = 2 resonance 
with total spin, parity J p =0~, isospin T=even, a resonance energy of M41— 2065 MeV, and 
a free decay width of T n NN m 0.5 MeV works extremely well in describing all available 
DCX data P,|10|. To exclude a possible nuclear structure explanation JTT|] , experimentalists 
have searched for the d' in proton-proton collisions, e.g. in p p — > p p tt + tt~ . A 4<r 
enhancement over the background exactly at the d' position has been observed in the ir~pp 
invariant mass spectrum (T2| . 



Two alternative models for the structure of the d! have been discussed in the literature. 
If the d! is predominantly composed of two colorless three-quark clusters, it could be a 
resonance in the NN* channel, where N* is the first negative parity excitation of the nucleon 
at 1535 MeV. Because the empirical mass of the d! is far below the N*N threshold, this 
requires an enormous binding energy for the bound NN* system. This interpretation is 
reminiscent of the deep- lying bound states in the Moscow NN potential ||13|| . The cf-dibaryon 



has also been interpreted as an isospin T = 2 NA resonance [fbl, a %NN state with T = 2 



TBI , and most recently as a 7rNN state in the isospin T = channel |L6|, located at 2018 
MeV. 

It is important to note that due to its quantum numbers (L = 1,S' = 1,T = 0), 
the d' cannot be composed of two ground state nucleons, because such a state is Pauli- 
forbidden. The Pauli principle demands that the quantum numbers of a two-nucleon system 
satisfy (— l) L+5+T = —1. Therefore, the d! cannot simply decay into two nucleons. A decay 
into two nucleons is possible only if a pion or a photon is emitted simultaneously. The 
correspondingly small phase space naturally explains the small width of the d! . 

On the other hand, quantum chromodynamics (QCD) does not preclude the possibility 
that the quarks in the initial ttNN system rearrange themselves into energetically favorable 
but unobservable colored clusters [JDJ. This idea is not new. Several scenarios for the 
arrangement of quarks and antiquarks into colored clusters combining to an overall color 
neutral B = 2 system have been proposed [ T8— pT|] . An early suggestion is the 'demon' 



deuteron, which is a J p = , T = dibaryon consisting of three pairs of diquarks [|22 



Other partitionings have been studied in the stringlike bag model |P3 , p4|| . This model 
predicts that a rotating (L = 1) dumbbell-like configuration with a colored diquark (g 2 ) and 
tetraquark (g 4 ) cluster connected by a flux tube of color-electric field lines is the state with 
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lowest energy for a six-quark system with d! quantum numbers. Fig. [I] shows how such a 
q 2 — q 4 cluster system could be formed from the initial 7rNN state. 

A shortcoming of the stringlike bag model p3|,|24|| is that it employs rigid impenetrable 



tetraquark and diquark clusters at the ends of a rotating color string. This model neither 
allows that the clusters merge into a single compound six-quark bag nor that they exchange 
quarks. Only the quarks within the individual clusters are antisymmetrized but not the 
quarks belonging to different clusters. Therefore, the stringlike bag model does not fully 
respect the Pauli principle. The stringlike bag model is a good approximation only for high 
angular momentum states (L > 5) In this case, the system is fairly elongated because 
of centrifugal forces, and the probability of cluster overlap is small. However, for a low lying 
L = 1 excitation, such as the d', one expects considerable overlap between the clusters and a 
substantial amount of quark exchange between them. We are sceptical that an L = 1 state 
is sufficiently stretched for the stringlike bag model treatment to be valid. The condition 
of validity is I > 2Rq PBfl , where / is the length and _R the radius of the color flux tube 
connecting the colored quark clusters. For the d', this condition is barely satisfied, and one 
anticipates corrections to the bag model prediction M# ~ 2100 MeV p3p4f . 

The purpose of this work is to study the dynamics of diquark-tetraquark relative motion 
using a microscopic quark model Hamiltonian and the Resonating Group Method (RGM). 
The model used here accurately reproduces the mass of the deuteron, which is the only 
established dibaryon. By comparing the RGM solutions with our previous results |26|-|28|, 
employing a single s 5 p 1 six-quark "bag" , we can test the validity of the assumption under- 
lying the stringlike bag model, namely that the d! is a stretched q 2 — q A system. In the 
present work, all complications arising from the quark exchange interactions (Pauli princi- 
ple) are rigorously taken into account, and their effect on the d! mass and wave function is 
investigated in some detail. In addition, the present paper critically examines the validity 
of the assumption of a universal two-body confinement strength and explores how a possible 
reduction of the effective confinement strength in a compound six-quark system affects the 
mass and the size of the d! . 

The paper is organized as follows. Section |T| provides a short description of the chiral 
quark model used. Section [Til] presents the six-quark Resonating Group Method approach. 
Section fV] discusses the numerical results for the mass and size of the d! for several confine- 
ment models and compares the present predictions to those of the single-bag shell model 
and the stringlike bag model. Our summary and conclusions are given in Section [V|. An 
appendix contains explicit expressions of the norm and Hamiltonian kernels needed for the 
solution of the RGM equation of motion. 

II. THE CHIRAL CONSTITUENT QUARK MODEL 

In the chiral constituent quark model a system of n-quarks with equal masses m q = 313 
MeV=mjv/3 (SUp(2)) is described by the Hamiltonian 

n 2 p2 n n 

H = ES + |9 - + E (r 4) r,) + E r,), (1) 
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where r;, are the spatial and momentum coordinates of the i-th quark and P is the total 
momentum of the n-quark system. The exact removal of the center of mass kinetic energy 
by the third term represents an important advantage of the present approach. 

The residual interactions V res = V OGEP + y OPEP + yOSEP k e t weeri quarks shown in Fig. 
model the most relevant properties of QCD, such as asymptotic freedom at high, and 
spontaneous chiral symmetry breaking at low energies. The one-gluon exchange potential 



V««(r„ r,) = ^A, ■ A,{i - i (l + 2 -a t ■ a,) f (r) j (2) 

provides an effective quark-quark interaction that has the spin-color structure of QCD at 
short distances. Here, r = — r^, and Oi is the usual Pauli spin matrix. The Aj represent 
the Gell-Mann matrices of SU(3) co i or . We neglect tensor and spin-orbit forces, which have 
been shown to be of minor importance for the d' ||27|j . Only central potentials are considered 
in this work. 

The spontaneous breaking of chiral symmetry of low-energy QCD by the physical vacuum 
is responsible for the constituent quark mass generation []30|], as well as for the appearence of 
pseudoscalar and scalar collective excitations of the vacuum (tt and o fields). These collective 
degrees of freedom couple directly to the constituent quarks. In the present quark model 
with two quark flavors, this mechanism is modelled by regularized one-pion and one-sigma 
exchange potentials between constituent quarks H3lH33|: 



fl 2 1 A 2 / p-m^r p -A n r\ 

vOPEP ^ *> - 1 12^1 Ai^r- ■ « ■ **< V - — > < 3 > 

q 7T 7T \ / 

V°^(r t ,r 3 ) = -^-rA^ {- — > W 

47r Az — mi \ r r 



% = m 2 ^(2m,) 2 + m 2 , = = A. (5) 

The nq coupling constant g 2 9 /47r is related to the 11N coupling constant g1 N /An by g wq = 
(3/5)(m q /mN)giTN with = 13.19 @]. The cut-off mass A describes the finite size of 



with 



the constituent quark due to its pion cloud 
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</ X 2 ' 



7* 2 — -777 . (6) 



Here, we choose A = 4.2 fm 1 which gives r q = 0.41 fm, and leads to a soft nN form factor 



35|j . The sigma parameters are fixed by Eq.(|^). In particular we use m a = 626 MeV. 

The usefulness and phenomenological success of the model has previously been demon- 
strated. For example, the chiral constituent quark model has recently been used to calculate 
the positive parity spectrum of the nucleon and various electromagnetic observables of the 



ground state octet and decuplet baryons and their excited states |35j . It has also been applied 
to calculate the properties of the deuteron |33]] , NN phase shifts [p2| , p3| , hyperon-hyperon 
interactions, as well as the H-particle [BHj. 
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A. Confinement models 



There are mainly two types of confinement potentials discussed in the literature, namely 
the two-body confinement potential introduced by Lipkin [37j], and the color flux tube 



(string) model of confinement. In the color flux tube model, the color string describes 
the hidden gluon degrees of freedom that are necessary to preserve color gauge invariance 
of the underlying field theory. In a many quark system, or in the interaction region of 
hadron-hadron collisions, the strings can change their positions and oscillate between dif- 
ferent configurations. This continuous flipping of color strings can be effectively described 



in the flip-flop model fl38fl , which has been introduced to avoid the long-range color van der 
Waals forces in hadron-hadron interactions. In the flip-flop model, the confinement interac- 
tion between any pair of quarks depends on the position of the remaining quarks; it thus 
contains many-body operators. 

Here, for reasons of simplicity, we consider several two-body confinement potential models 
of Lipkin-type, which differ in their radial dependence (see Fig. 3) but not in their color 
structure. We should keep in mind that these potentials may not be adequate for describing 
the complicated dynamics of changing color strings in a many-quark system. 

In the constituent quark model one often takes a simple quadratic confinement potential: 

V; conf (r 4 ,r J ) = -a^A J -A J (r 4 -r J ) 2 . (7) 

Later we use this model to argue that the usual assumption of a universal two-body con- 
finement strength a® is too restrictive and may not be adequate for compact six-quark 
systems. 

Lattice-QCD calculations find that the quark-antiquark potential of QCD is linear: 

Vr i (r l ,v J ) = -a^X t .\ 3 \(r l -r J )\, (8) 
The linear form has also been used in our previous calculation p6|p7| using the Translation- 



ally Invariant Shell Model (TISM). 

Another interesting parametrization is the r 2//3 confinement potential: 

K conf (r J ,r,) = -aWA,-A J |(r l -r J )r /3 . (9) 
In Ref. [p9| it is shown that a r 2 / 3 confinement potential and a color-Coulomb potential leads 



in the framework of the nonrelativistic Schrodinger equation to the observed linear Regge 
trajectories of hadron masses. 



Finally, we consider the color-screened error-function confinement [40 

2 rx 



i^Xi ■ A,erf(//r), erf(x) = -= / e~ z ' 'dz . (10) 

\/7r Jo 



Vr i (r l ,r J ) = 

This potential rises linearly for small r, but as a result of quark-antiquark pair creation 
(color screening) grows only weakly for intermediate r and finally goes to a constant value 
at large r (see Fig. 3). Lattice calculations show [^] that such a behavior of the effective 
quark-quark potential arises if quark-antiquark loops are taken into account. The inverse of 
/i is called color- screening length for which we take = 0.8 fm. This potential has recently 
been used by Zhang et al. fl41j . Alternatively, an exponential form V^°^ = (1 — exp(— /xr 2 )) 



that matches the r 2 behavior for small r has been suggested and used by Wang et al. [02] 
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B. Determination of the parameters 



For the nucleon and A wave functions, we use s 3 harmonic oscillator ground state wave 
functions $at(a) 

/ i \ 3/2 

I **(A) >= exp (- (p*/(4bl) + X 2 /(3bl))) | ST > N ^ x | [111] (11) 

where the Jacobi coordinates p and A are defined as p = r± — r 2 and A = r 3 — (ri + r 2 )/2. 
The harmonic oscillator parameter 63 describes the mean square (ms) matter radius of the 
nucleon. The matter radius is defined with respect to the center of mass coordinate of 
the nucleon and is given as 



6 i=i 



$ TV = b 2 3 . 



As usual, the parameters of the Hamiltonian of Eq. ([!]) are fitted to the nucleon and A 
ground state masses. The masses of the N and A are calculated as expectation values of 
the Hamiltonian of Eq. ([I]) between baryon ground state wave functions $jv(A)- Imposing 
=313 MeV and A=4.2 fin -1 , the parameters af\ a s , and 63 are determined by the 



m 



requirement that the iV(939) and A(1232) masses are reproduced, and that the nucleon 
mass is stable with respect to variations in 63: 

M N (b 3 ) = 3m„ = 939 MeV, M A - M N = 293 MeV, dM ^ h ^ = . (12) 

Although the three coupled equations (|P2| ) for the parameters af \ a s , and 63 are non- 
linear, a^> is predominantly determined by the requirement that the nucleon has its physical 
mass Mtv = 939 MeV, a s by the mass difference Ma — Mjv, and 63 by the stability condition 
dMN(bs)/db3 = 0. The parameters for the different confinement models as determined 
by Eq. ( |T2"D are shown in Table |I[ Parameter sets I-III are obtained with the quadratic 
confinement of Eq. ([FD and sets IV,V,VI are obtained with the linear, r 2//3 , and error-function 



confinement of Eqs.(^-TU) respectively. 

One can qualitatively understand the trend in the values of the parameters a s , and 
63 for the different confinement models by considering the sequence quadratic (set II), linear 
(set IV), r 2 / 3 (set V) and error-function (set VI) confinement for the case of the nucleon (see 
Fig. |H). One observes that the stronger the confinement potential grows at small distances, 
the further out the nucleon wave function (larger 63) extends, in order to minimize the 
nucleon mass. With increasing 63, the importance of the one-pion-exchange interaction 
decreases, and the contribution of the one-gluon-exchange to the NA mass splitting, -hence 
a s - must increase correspondingly. Finally, because the kinetic energy contribution to the 
nucleon mass is lower for larger nucleon sizes (63), a bigger confinement strength is 
needed to fit the experimental nucleon mass Mtv=939 MeV. 
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III. A COLORED DIQUARK-TETRAQUARK CLUSTER MODEL OF THE D' 

DIBARYON 



In the DCX reaction, the incident pion and the two correlated nucleons may come suf- 
ficiently close for the quarks to make use of the additional possibilities provided by their 
color degree of freedom. In the interaction region, there are several possibilities to make an 
overall color-singlet. The quarks could clusterize into two color singlets, two color triplets, 
two color sextets, and two color octets: 

(lOl) 1 ; (3&3) 1 ; (6 6) 1 ; (8 <g> 8) 1 . 

We recall that a q 3 — g 3 clusterization of the d! with colorless three-quark clusters is either 
not allowed by quantum numbers (Pauli principle) or by energy considerations. For example, 
a iV(939)iV(939) clusterization with relative angular momentum L = 1 is not allowed by the 
Pauli principle. Similarly, an NA system with L = 1 can only have isospin T = 2, while 
present experimental results seem to favor the T = assignement. Furthermore, an L = 
iV(939)iV*(1535) system, where the negative parity resides inside one cluster, is far above 
the experimental d! mass. A recent quark cluster model calculation has found no evidence 
for a deeply bound state in this system pE3]| . Therefore, a iV(939)iV*(1535) dibaryon decays 
strongly into NNit. It has a large width and thus cannot be a viable candidate for the df. 
Thus, the d' is presumably not a state which is predominantly composed of two colorless 
three-quark clusters. 

In the stringlike bag model, the color-triplet (3) tetraquark, color-antitriplet (3) diquark 
clusterization is the energetically most favorable configuration for a system with d! quantum 
numbers. For such a system, a relatively small mass of ~ 2100 MeV has been obtained 



The six-quark wave function for the cf -dibaryon is expanded into the 3 ® 3 tetraquark- 
diquark cluster basis 



$? T - 1,TT -V.Ar,»7T) x [2H]? 



J=0,T=0 

C=0 



[222]^ , (13) 



where $^ T_1 ' TT_0 (p T , \t,t) t ) and $^ D_0 ' Ti3_c \p D ) are the internal wave functions of the 
tetraquark (T) and diquark (D) clusters, respectively, and xl=i(JI) is the yet unknown 
relative wave function of the two colored clusters, projected onto good angular momentum 
L=l. The Jacobi coordinates introduced in Eq. are depicted in Fig. |[ As usual 
in Resonating Group Method (RGM) [45] calculations, ground state s 2 and s 4 harmonic 
oscillator wave functions are used for the diquark and tetraquark clusters, respectively. 

Eq. ( |T3"D shows that the color triplet (3) tetraquark with mixed permutational symmetry 
[211], and the color antitriplet (3) diquark with permutational symmetry [11] is coupled to 
an overall color-singlet (1) six-quark state with mixed permutational symmetry [222] in color 
space. In contrast to mesons and baryons, the color state of a genuine q 6 system, is by itself 
not fully antisymmetric. There are both antisymmetric and symmetric quark pairs. Thus, in 
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a compact g 6 system, there is no factorization of the color space and the flavor-spin-orbital 
space, i.e., the color dynamics and the orbital structure of the system are correlated in a 
more complex manner ||37|| . This can be seen as an indication that the transition between 
q 3 and compact q 6 systems is not as trivial as usually assumed. 

We recall that the stringlike bag model employs impenetrable diquark and tetraquark 
clusters, and the clusters cannot merge into a single bag. In contrast, an RGM calculation 
allows for a continuous transition from the compound q 6 state, where all quarks are in a single 
potential well, to the q 2 — q 4 clusterized state, where one has two clearly separated bags. 
There is no artificial boundary between these extreme configurations, and they are both 
described by one and the same Resonating Group Method wave function. Furthermore, the 
q 3 —q 3 and q^—q 1 partitioning into colored clusters, as well as the q 3 —q 3 split into color-singlet 
clusters is automatically included in the present theory. These important properties are a 
consequence of the Pauli principle on the quark level which is ensured by the antisymmetrizer 
A 

A=l- 8P° STC + QPg STC Pg STC , (14) 

where P° STC is the permutation operator of the i-th and j-th quark in orbital (O), spin- 
isospin (ST) and color space (C). The direct, one-quark- and two-quark-exchange contribu- 
tions associated with the three terms in the antisymmetrizer of Eq. (|14"D are shown for the 
case of the one-gluon-exchange potential of Eq. (|2|) in Fig. ||. 

The solution for the unknown relative wave function Xl(R-) an d the unknown eigenenergy 
E = Md' is obtained from the Ritz variational principle 
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(V d ,\H-E\V d ,) 
(*d'|*d'> 



0, (15) 



where the variation is performed with respect to the relative wave function x^(R). After 
expanding the relative motion wave function Xl(R-) i n a finite sum of Gaussians centered 
around the generator coordinates Sj (Generator Coordinate Method (GCM)) 

/ 4 \3/4 / o \ 
Xl(R) = Pl £a (^f) exp^-iR-stf/bl) =► I**) =:AY,C t |GCM,(s,)) 

(16) 

and projecting onto good angular momentum L via Pl = l/v47r / dsi ■ Y L (§i) (cf. the ap- 
pendix for more details), Eq. (p~5| ) transforms into a generalized algebraic eigenvalue problem 

Y^n^sj) c { = e Y^iM,^) Ci . (17) 

i i 

The analytic expressions for the Hamiltonian kernel TCij(si,Sj) = (GCMjl^Ai/IGCMj) and 
the norm kernel Mij(si,Sj) = (GCMj|^4|GCMj) of Eq. ( |TTD can be found in the appendix. 
The norm matrix J\fij(si,Sj) on the r.h.s. of Eq. (|i~7D reflects the non-orthogonality of the 
Gaussian basis functions in Eq. (0). For simplicity, and to avoid a proliferation of param- 
eters, the same variational parameter b§ for the internal and relative motion wave functions 
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of the six-quark system is used. This greatly facilitates the calculation of the norm and 
Hamiltonian integrals and restricts the variational space spanned by Eq. ([if]) only slightly. 

Having solved the generalized eigenvalue problem, we determine the mass of the d! and 
the optimal harmonic oscillator parameter 6 6 , which plays the role of a nonlinear variational 
parameter, by minimization of the d' mass: 

dE d ,(be)/dh = 0. (18) 

We emphasize that the b$ values obtained in this way are some 20-30% larger than the 
corresponding values 63 which minimize the nucleon mass. This observation is the basis 
for our conjecture that the two-body confinement strength in a compound six-quark system 
such as the d! is weaker than in a three-quark system (see Sect. |1V C|) . 



IV. RESULTS AND DISCUSSION 



In this section we present the results of an RGM calculation for the mass and size of 
the d! dibaryon and compare it with the corresponding shell model and stretched bag model 
results HU. Because a corresponding RGM calculation for the deuteron, based on the 
same Hamiltonian and colorless nucleon cluster wave functions accurately reproduces the 
deuteron energy and wave function |14| , we are confident that the absolute mass scale of our 
predictions is correct, and the calculated d! masses are reliable. 



A. Mass of the d' 

The mass of the d! is obtained as the lowest eigenvalue E of the bound state GCM 
equation (|17D, and by subsequent optimization with respect to b 6 according to Eq. fll8|) . It 
is interesting to isolate the effect of the Pauli principle on the mass of the d! by performing 
a calculation where the quark exchange kernels have been turned off, i.e. only the first term 
in the antisymmetrizer of Eq.(|l4]) is taken into account, before solving Eq . (|17|) . We use 
the same b$ as in the full calculation including all exchange kernels. Results for and 
b§ are shown in Table || for the parameter sets of Table [[[ A comparison of the results for 
the d! mass with and without quark exchange show that the quark exchange interactions 
contribute an additional energy of 50-120 MeV depending on the model of confinement (see 
sets I- VI in Table |TJ). 

We note that the d' mass calculated with model I (quadratic confinement and one-gluon 
exchange) is some 550 MeV above the experimental mass of the d! . An RGM calculation 
without 7r- and cr-meson exchange forces in the Hamiltonian, but with a q 3 -q 3 clusterization 



of the quarks has already been performed many years ago j|7fl • The quark-quark interaction 



of Ref. [|7[ is very similar to model I, and the method of fitting the quark model parameters 



to single-baryon properties is the same as in the present work. Although the authors of Ref. 
p7| use two color-octet three-quark clusters, i.e., (8C5 8) 1 whereas we use color-(anti)triplet 
diquark and tetraquark clusters, i.e., (303) 1 , it is nevertheless meaningful to compare their 
and our results. 
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Due to the quark antisymmetrizer A of Eq. ([TJ|) acting on all six quarks, the three 
different arrangements of the six-quarks into colored clusters, namely q 3 -q 3 , q 2 -q 4 , and q x -q 5 
are approximately equivalent for the relevant intercluster distances involved. In fact, in the 
region of complete cluster overlap, that is, in the limit for s, = Sj — > all three partitions 
have the same six-quark shell model limit, i.e. the unique s^p 1 TISM state. Thus, by virtue 
of the antisymmetrizer, the RGM wave function not only contains the (3 (g) 3) 1 component, 
but also the (ScgiS) 1 component. It is reassuring that the results of Ref. H7| Mrf/=2540 MeV, 
and of the present work M^/=2610 MeV (see set I in Table [H|) are in very good agreement. 
However, both are ~ 500 MeV above the experimental resonance position. Before drawing 
any conclusions concerning the existence of the d', we first study how these results depend 
on the choice of the quark-quark interaction. 

In model II, we use the full Hamiltonian including one-gluon exchange, tt- and a-exchange 
BTjfl , and a quadratic confinement potential between the constituent quarks. The inclusion 



of the chiral tt- and cr-exchange interactions reduces the d' mass by 160 MeV compared to 
model I. It is, however, with 2440 MeV still about 380 MeV higher than the experimental 
value. 

Model III is very different from all other confinement models. It uses a different confine- 
ment strength for the three-quark (nucleon) and the six-quark system (d f ) and leads to a d' 
mass close to the experimental result. Before we try to justify this apparently 'adhoc' step, 
we study the dependence of the d' mass on the radial form of the confinement potentials 
given by Eqs. - [TOD and shown in Fig. 3. 

The slope of models (IV-V) is for r > 1 fm considerably flatter than that of the quadratic 
confinement of Eq. ([?[). For these potentials, the stability condition Eq.(^) leads to larger 
oscillator parameters b§ (larger dibaryons), which in turn results in a reduction of the kinetic 
energy and the d' mass. For example, the linear confinement potential (model IV) reduces 
the d! mass by about 90 MeV in comparison to the quadratic confinement (model II). For 
the r 2 / 3 confinement (model V), the corresponding reduction is 130 MeV. 

Model VI, the color-screened error-function confinement, displays a qualitatively different 
behavior at small and large interquark distances. At short distances, it is approximately 
linear and at larger distances it goes to a constant. While baryons feel mainly the linear 
rising part, dibaryons are due to their larger size also very sensitive to the long-range tail 
of the confinement potential. The screening of color charges at larger interquark distances 
leads to a reduction of the d! mass by 150 MeV compared to the quadratic confinement 
(model II) and also to a larger value for the oscillator parameter 6 6 .Q However, the d! mass 
is with 2288 MeV still some 220 MeV above the experimental value of M^ p = 2065 MeV. 
If we use the 1 — exp(— fir 2 ) color-screening model of Ref. with /i = 1 fm -2 for all quark 
pairs in the system, we obtain Ma> = 2468 MeV at = 0.82 fm. In this context, we mention 
that at interquark distances corresponding to 3-4 times the screening length l//x= 0.8 fm, 
the effective confinement force — dV^ (r)/dr between the colored clusters tends quickly to 
zero. The stability of the system is nevertheless preserved due to Eq . fli~8|) and our use of 



x In the suggested range for the color-screening length 0.8fm < 1/fi < 1.2frn [41|, our results are 
not very sensitive to the actual value of 
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Gaussian trial wave functions. Model VI represents a limiting confinement model especially 
for dibaryons, and we expect that the d' mass cannot be further reduced by an significant 
amount by considering another radial form of the confinement potential. 

In Table [H] we list the individual contributions of the kinetic, confinement, one-gluon-, 
one-pion-, and one-sigma-exchange potentials to the d! mass for the different confinement 
models. The results in parentheses give the corresponding contributions without the quark 
exchange diagrams. One observes a drastic reduction of the kinetic energy for the larger 
dibaryon systems. We also see that the dominant attraction is due to the color-Coulomb 
part of the interaction. 

Thus, we observe a clear trend when going from a quadratic (model II) to a color-screened 
confinement model (VI): The slower the increase of the confinement potential at larger r, the 
larger the size 66, and the smaller the mass of the d! . However, even for the limiting case of 
the color-screened confinement potential, the results for the d' mass are 220 MeV above the 
empirical d' mass. We conclude that in the present constituent quark model framework, a 
change in the radial form of the confinement potential is insufficent to generate a cf-dibaryon 
with a mass of 2065 MeV. 



B. Wave function and size of the d' 



We turn now to the discussion of the relative wave function as obtained from eigenvalue 
equation (|T7|) and the radius of the d' . Fig. |^ shows the wave function xl=i(R) for the two 
extreme confinement models II and VI, using the same Hamiltonian as in the three-quark 
sector. The RGM wave function is calculated both with (Pauli-on) and without (Pauli-off) 
the quark exchange diagrams. For the quadratic confinement model II, one observes that the 
relative wave function has already died out at intercluster distances of about 2.5 fm, while 
for confinement model VI it still has an appreciable amplitude even at 4 fm. The quark 
exchange diagrams lead for model II to an additional attraction between the two colored 
clusters, as reflected by the smaller radial extension of the relative wave function, while they 
are rather unimportant for model VI. 

A quantitative measure of the extension of the relative cluster wave function is the root 
mean square (rms) distance, R§ GM , which measures the mean distance between the q 2 and 
q A clusters 



/ r>RGM\2 
[ K d> ) 



(xl=i(R)\R 2 \xl=i(R)) 
(xl=i(R)\xl=i(R)) 



(19) 



Here, R is the relative coordinate connecting the centers of the two clusters. This distance 
corresponds to the length I of the color flux tube in the stringlike bag model. 

The rms distance of Eq. ( |I~9"D does not yet take into account the finite size of the clusters. 
In the present model these are given as 



r 2 



D 



i=l 



R 2 



D 



4 



T 



4f 



Rr) 2 



i=l 
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b 2 . (20) 



Here, R^ and R T are the center of mass coordinates of the di quark and tetraquark, respec- 
tively. One then finds for the total d! radius in the RGM formalism 
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r 2 D + r> T + (~Rr M f + r 2 q , (21) 



where the last term is due to the finite size of the constituent quark. In Table fV|, the radii 
of the diquark, tetraquark, and the relative cluster wave function, R^ GM (with are shown 
for different confinement potentials. Neglecting quark exchange between clusters increases 
R$ GAI typically by about 10% (numbers in parantheses) . This shows that the quark exchange 
matrix elements provide additional color attraction. Depending on the confinement model, 
the rms distances R^ GM lie between 1.0 and 2.0 fm, while the corresponding tetraquark 
radii, r?, vary between 0.74 fm and 1.15 fm. 

In order to have another measure of the size of the d', we calculate the rms distance 
between an arbitrary pair of quarks in the dibaryonQ 



C) 2 = ** 



^Vpairs i < j 



* d' ) (22) 



where N pairs = 15 is the number of quark pairs in a six-quark system. The spatial extension 
of the quark distribution in the d' is close to the value given by Eq. fl2"2]) . The radius defined 
in Eq.(|22|) facilitates the comparison with the shell model description of the d! []26| , |27 , 46 . 
The numerical results for r^ GM are given in the second to last column of Table |TV| . 

1. Comparison with the six-quark shell model 

In Fig. |6| we also compare the q 2 — q A relative RGM wave function with an unperturbed 
six-quark lhu (s 5 p 1 ) harmonic oscillator wave function 

Xl=i(R) = i = R exp --T2- , (23) 



with the same oscillator parameter 6g. Our previous results |27| using the translationally 
invariant shell model (TISM) have shown that admixtures of excited shell model states to 
the unique lowest lying s 5 p 1 , lhui shell model state of Eq. (p3|) are small. 

For the quadratic confinement model II, the d! wave function does not display a pro- 
nounced clusterization but resembles rather closely the six-quark harmonic oscillator wave 
function of Eq . (p3|) . The agreement between the pure harmonic oscillator wave function and 
the RGM wave function is more complete when the quark exchange diagrams are included. 
On the other hand, the RGM wave function of model VI extends to larger intercluster dis- 
tances R than the unperturbed shell model state of Eq.(^). This is a consequence of the 
color-screened confinement potential, which provides only weak confinement forces at large 
interquark distances. In order to describe the long-range RGM wave function in the shell 
model, one would need substantial admixtures of excited harmonic oscillator states. As one 



2 Note that in the case of the nucleon the mean distance between any pair of quarks is 63 • y/3 — lfm. 
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might expect, the effect of the quark exchange diagrams on the wave function is here much 
smaller than for the quadratic confinement case, which can also be seen from Table ||. 

For the wave function of Eq. (^), the mean square (ms) distance between the clusters 
is given by 



{RTf = {xt\ R 2 x££) = f K (24) 

For all but one model, the rms distance between the q 2 and q 4 clusters as calculated in RGM is 
slighly larger than the corresponding quantity in the harmonic oscillator model. Only in the 
case of the color-screened confinement potential (model VI), we find a pronounced increase 
of R$ GM by 40% compared to R$°, indicating that the system is partially clusterized. 
In the case of the TISM wave function including shell model configurations up to 3tku 
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27J, the rms-distance between any two quarks in the d! is given by 



r 



I ISM ^2 _ ; 2 / i / 1 



qq , -o V 5 ■ v- - / g 



h\ - + (l-a\ , (25) 



where a 2 denotes the probability of finding the (unique) energetically lowest lying N=l shell 
model configuration in the total TISM wave function of the d! . The values given in Table 



IV] are obtained for a harmonic oscillator parameter bQ, which minimizes the d! mass (see 
Table |IJ. 

For models I - V, the rms distance between an arbitray pair of quarks r^ M follows 
closely the corresponding TISM result r^ SM . Only for the color-screened confinement (set 
VI), the average distance between quarks in the cluster model r^. GM is increased by about 
10% with respect to the corresponding shell model result . This is mainly due to 

the intercluster quark pairs, which due to the screening are preferably at larger relative 
distances. 

Thus, a comparison of the cluster model and shell model results |27| reveals an overall 
agreement for the orbital structure and size of the d! . In addition, the two calculations 
agree also for the individual kinetic and different potential energy contributions to the d! 
mass (as given for the cluster model in Table [Hl|) , and thus also for the d! mass itself. The 
agreement is quantitative for models I-V. However, the longer tail of the RGM wave function 



of model VI shown in Fig. |6| and the larger rms-radius listed in Table [IV| show that certain 
confinement models lead to deviations from a pure s^p 1 shell model structure and favor a 
clusterization of the quarks into colored quark clusters. 

Let us briefly discuss the reasons for the quantitative agreement between the two struc- 
turally quite different calculations. The outer product of the [4] G (tetraquark) and [2] G 
(diquark) orbital symmetries (s 4 and s 2 ) gives according to Littlewood's theorem the fol- 
lowing Sq permutational symmetries in orbital space 

[4] ® [2} = [42] © [51] © [6] . (26) 

For d! quantum numbers, the fully symmetric [6]o orbital symmetry, corresponds to a spu- 
rious center of mass excitation of the six-quark state, and is automatically excluded in both 
approaches. The orbital symmetries [51]o and [42] o in Eq.(refsym) are also included in the 
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enlarged N=3huj shell model basis p7| . Analogously, the outer product of the two clusters 
in spin-isospin space leads to 



[31] sr ® [2]st = [51]st © [42] sr © [33] sr © [411] sr 
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[321] ST . (27) 
shows that the q 2 — q 4 



Comparison with Eq. (10) in our previous TISM calculation 
cluster model wave function comprises the same S^-symmetries in spin-isospin space (with 
the exception of the [221 1}st symmetry which does not appear in the cluster model) as 
the enlarged shell model basis. Thus the trial function space spanned by both sets of basis 
functions is nearly equivalent. 



2. Comparison with the stringlike bag model 

It is interesting to compare the description of the d' in the cluster model with the one in 
the stringlike bag model. With respect to the sizes, the cluster model values R§ GM and vt 
have to be compared to the length I and radius Rq of the color flux tube in the stringlike 



bag model. In the stretched bag model |£J , the radius R and length I of the flux tube are 
given as 

respectively, where a — 1.1 GeV 2 is the "universal" string tension, B = 59 MeV fm -3 the 
bag constant, and = 2100 MeV is the mass of the resonance. One then obtains for the 
"universal" radius of the color flux tube Rq = 1.1 fm and I = 2.3 fm |[23||. 



According to Ref. ||23|| , the stretched bag model is only valid if / > 2Rq, and for a low-lying 
resonance of Mj> = 2100 MeV this condition is barely satisfied. If we now identify the bag 
model quantities I and i? with the cluster model quantities R$ GM and r T respectively, we 
see that this condition is generally not satisfied in the present calculation. In fact, from table 



rV| it is evident, that the sum of the cluster radii exceeds the intercluster radius by about 
half the diquark radius. This means that there is considerable overlap between the clusters. 
Only for the color-screened confinement potential, we observe a modest clusterization that 
is, however, insufficient to justify a stringlike bag model treatment. 

At this point, we comment on the universality of the string tension a in the stringlike 
bag model ||23|| , which predicts a series of orbitally excited meson, baryon, and dibaryon 
states with orbital angular momentum L, and with masses given by 

M 2 = aL + M 2 . (28) 

The mass of the clusters at the ends of the string, M , is calculated in the spherical bag 
model, and a is the common string-tension for mesons (q — q), baryons (q 1 — q 2 ), and 
dibaryons (g 4 — q 2 ) which is given by: 



a 



SrrBaJ 2 , (29) 



where B is the bag constant determining the constant color-electric field strength E in the 
string, a s is the quark-gluon coupling, and is the eigenvalue of the quadratic Casimir 
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operator of SU(3) co i or in a given color representation. It is important to recall that the 
string tension a and the related transverse radius of the string R are defined only for well 
separated clusters with sufficiently high relative angular momentum L [fj3| . 

A universal string tension a results because it is assumed that the individual clusters are 
so far apart that they are always in the color-triplet 3 or color-antitriplet 3 representation, 
for which the Casimir operator has the common eigenvalue = 4/3. In this case, there is 
only one type of color string with a unique string tension given by Eq. . Obviously, this 
assumption is justified as long as the clusters do not overlap. However, we have seen that 
for most confinement models there is considerable overlap between the colored diquark and 
tetraquark clusters. Consequently, the tunneling of quarks from one cluster to the other 
can change the color-representation (e.g. into 8 <8> 8, or 6 ® 6) of the clusters and therefore 
also the string tension. For the octet and sextet color representations, which are not present 
in mesons and baryons, one has = 3 and f@ = 10/3 respectively. Thus, one obtains a 
stronger string tension in both cases. The authors of Ref. [^] only include the energetically 
most favorable case of a 3 ® 3 string. If the Pauli principle is to be rigorously satisfied in 
Ref. |23[| also the 6 ® 6 and the 8 (g> 8 color configurations would have to be admixed and the 
corresponding d! mass would come out higher. This qualitatively agrees with our results for 
the d' mass (see Tables |T|-|TT| ) which are increased by the quark exchange interactions. 

There is another difference between the present RGM calculation and the stringlike bag 
model. For the meson (q — q), baryon {q — q 2 ), and dibaryon [q A — q 2 ) resonances, the sizes of 
the color sources (clusters) and therefore the transverse radii of the flux tubes are in reality 
rather different. A tetraquark in a dibaryon has a larger radius than a single quark in a 
baryon. Although a tetraquark and a single quark carry the same color charge, the energy 
density of the color electric field, E 2 /2, in the neighborhood of an extended tetraquark is 
smaller than in the neighborhood of a single quark. This difference is rather important for 
low angular momentum states. In the stretched bag model, the energy density of the string 
is given by the bag constant B. A decrease of the energy density of the string for g 4 — q 2 
dibaryons compared to q — q 2 baryons is related to a decrease of the effective bag constant B 
f\. According to Eq.(p9|) this leads to a reduced effective string tension for dibaryons. This 
reduction of the string tension clearly outweighs its increase due to quark exchange. 

In conclusion, the assumption of a universal string tension is not justified for ground 
states, such as the pion, nucleon, or the d' dibaryon. We emphasize once again that the 
stringlike bag model is valid only for higher angular momentum states, where the clusters 
are well separated due to the centrifugal barrier, and where the tunneling of quarks, as well 
as the size difference of the color sources can be neglected. 



3 The universality of the bag constant for baryons and dibaryons has been questioned before B 
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C. The d! as an indication for a weaker confinement strength in a compound 

six-quark system? 



So far we have studied how different confinement models affect the mass and characteris- 
tic size of the dibaryon. In all cases, we obtained dibaryon masses some 200-400 MeV above 
the experimental value. We point out that up to this point, we have used the same param- 



eters as previously determined from the properties of color-singlet baryons (see Eq.(12)). 
Although it is reasonable to expect such a universality of parameters for the microscopically 
better founded gluon-, pion- and sigma-exchange interactions, the assumption that the con- 
finement strength remains unmodified when going from a three-quark to a six-quark system 
is most likely too restrictive. 

At present, there is no commonly accepted theory of color confinement. Different models 
of confinement have a limited range of validity. From the experimental information on 
the excited nucleon spectrum, combined with certain assumptions concerning the Lorentz 
structure, radial form, and color-dependence, phenomenologically successful confinement 
models for baryons have been constructed. However, it is not at all obvious that one can 
extrapolate the experience gathered in the color-singlet three-quark sector to compound 
six-quark systems. We recall that the confinement strength a c is usually determined from 
the experimental level-spacing uj between excited single-baryon states. Unfortunately, there 
is hardly any three-star data on excited six-quark resonances that would shed light on the 
confinement dynamics in a six-quark system. 

We have seen that in the d', the average distance between any two quarks is somewhat 
larger than in the nucleon. Consequently, the d' tests the confinement interaction at larger 
quark-quark distances where new physical phenomena come into play. For example, at 



large interquark distances, quark-antiquark pair creation plj|42| leads to a screening of the 
original color charges. This color-screening effect can be simulated by a reduced effective 
confinement strength in the six-quark system as compared to the three-quark system. 

Although we do not know how to calculate the effective confinement strength for three- 
and compound six-quarks systems from first principles, we can gain some qualitative under- 
standing of the changing confinement dynamics using the harmonic oscillator model. If the 
effective quark-quark interaction were a pure harmonic oscillator confinement force, there 
would be an inverse proportionality between the confinement strength 

*° = ^fcv < 3 °> 

and the fourth power of the harmonic oscillator parameter brm- 

In model II, we have used one and the same two-body confinement strength for three- 
quark and six-quark systems in the numerical calculation, even though the variational prin- 
ciple of Eg. (fl^j and Eg. (18) tells us that the characteristic sizes of the d! (b e ) are about 



30% larger than those of the nucleon (63) . Thus, in a model where the two-body harmonic 
oscillator confinement is the only quark-quark interaction, there is a unique relation between 
the two-body harmonic confinement strength and the size of the system. We take this ob- 
servation as the basis for our conjecture that the effective two-body confinement strength 
could be weaker than in a three-quark system f4"6]| . Although the actual dependence of 
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the effective confinement strength on the size of the system may be somewhat different for 
more realistic confinement models, we expect that an inverse proportionality between the 
confinement strength and the size of the system remains. 

Model III differs from model II only in the strength of the parameter af£\ for which we 
have taken the value = 5.0 MeV/fm 2 « a\ 3 > /3 ||53|| . All other parameters in the Hamil- 
tonian are identical to the ones of model II. The proposed weaker two-body confinement in 
a six-quark system leads to a larger d! and a mass that is close to the experimental d' mass. 
Conversely, one could take the empirical d' mass as evidence for a weaker effective two-body 
confinement strength in a compound six-quark system because in the present model this is 
the only way to obtain the d' resonance mass, which is needed to fit the DCX data. 



D. Two-baryon vs. six-quark description of dibaryons 

The use of a reduced confinement strength in a compound six-quark system, where all 
quarks are basically in a single bag, is presumably not in conflict with our previous results 
for the deuteron. The deuteron is well decribed by an RGM wave function that is built from 
two colorless three-quark clusters. The contribution of compound six-quark states to the 
deuteron wave function is rather small. 

Here, we would like to contrast the description of the d' dibaryon as a compound six- 
quark state with the successful picture (in terms of quarks) of the deuteron as bound system 
of two color-singlet objects. This comparison clearly shows that the orbital structure of 
these systems is determined by an interplay of the Pauli principle and different parts of the 
effective quark-quark interaction. 

It has been known for more than a decade that quark exchange between nucleons (Pauli 
principle) together with the spin-dependent quark-quark interactions provide an effective 
short-range repulsion in the NN system [|53| . The latter prevents appreciable cluster overlap 



and is therefore mainly responsible for the successful description of the deuteron as a system 
of two colorless nucleons with an average distance of about 4 fm between the nucleons. To 
explain this in terms of a six-quark shell model we form the outer product of two s 3 nucleons. 
One then obtains the following spatial permutational symmetries of a six-quark system 

[3]o ® [3]o = [6]o © [42]o © [51]o © [33] G . 

In a six-quark system with even orbital angular momentum, only the [6]o and [42] o symme- 
tries occur. The important role of the orbital [42]o symmetry for the short-range repulsion 
in the NN system has been elucidated in Ref . |J3J . In the shell model, which is appropriate 
for short NN distances, the short-range repulsion is achieved through a destructive inter- 
ference of the excited s 4 p 2 six-quark state and the s 6 six-quark state. The spin- dependent 
quark-quark interactions lower the energy of the excited s A p 2 state and raise the energy of 
the s 6 state. As a consequence, the excited s 4 p 2 state is admixed with equal weight but 
opposite sign to the s 6 ground state so that it interferes destructively. This in turn leads 
to an almost complete cancellation of the wave function in the region of cluster overlap 
and consequently to the suppression of the deuteron wave function at short distances ||55| . 



Furthermore, explicit calculation shows that the deuteron wave function is insensitive to the 
details of the confinement mechanism. 
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On the other hand, in the case of the AA-decoupled d'-dibaryon, the radial form of the 
confinement potential crucially affects the mass and wave function of the system as we have 
demonstrated. The spin-dependent interactions are of minor importance. The latter point is 
reflected by the small admixtures of excited states to the lowest lying s 5 p x shell model state. 
Because the Pauli principle prevents a clusterization of the six-quarks into two colorless 
nucleons, the energetically lowest configuration is a compound six-quark state. Thus, for 
both, the deuteron and the d', the Pauli principle determines whether the residual spin- 
dependent or the confining interactions prevail, and hence determines the orbital structure 
of the system. 



V. SUMMARY 

In the present work, we have calculated the mass and wave function of a J p =0~, T=0 six- 
quark system, called d', in a colored diquark-tetraquark cluster model using the Resonating 
Group Method (RGM). This method determines the orbital configuration of the six-quark 
system dynamically, i.e. according to the given Hamiltonian. Thus we may test the validity 
of the assumption that the d! is a stretched diquark-tetraquark system. In contrast to the 
stringlike bag model, which employs a single non-antisymmetrized q 2 — q A dumbbell-like 
configuration, the present RGM calculation also includes other clusterizations, such as the 
q 1 — q 5 , q 3 — q 3 and the single q 6 state. This remarkable property is a consequence of the 
proper antisymmetrization of the total six-quark RGM wave function. 

A major purpose of this work has been to study the effect of quark exchange interactions 
between the colored clusters on the mass and wave function of the d' . The quark exchange 
interactions between the colored diquark and tetraquark clusters increase the d' mass by 
some 100 MeV and slightly decrease the size of the d' compared to a calculation without 
quark exchange. Our results for the d' mass and wave function are rather similar to a previous 



six-quark shell model calculation |2qj27|j46|] . This means that for confinement models I-V, the 
dl is in reality not a clusterized q 4 — q 2 state, but better described as a single six-quark s 5 p 1 
shell model state, where all quarks move in a common potential well. For this conclusion, 
it is crucial that the quark exchange diagrams are included in the cluster model. 

In order to investigate how our results depend on the model of confinement, we have 
studied various, commonly used confinement potentials. All models yield d' masses substan- 
tially larger than the empirical value. Depending on the confinement model we obtain d' 
masses that lie roughly 200-400 MeV above the experimentally required resonance energy of 
Md/=2065 MeV. Thus, we conclude that irrespective of the radial form of the confinement 
potential, we cannot describe a J p = 0~ T = dibaryon with a mass of 2065 MeV if we 
use the same confinement strength as in a three-quark system. Can one conclude that the 
d' does not exist? 

In a six-quark system, such as the d', the confinement interaction is tested in a heretofore 
unexplored regime where new physical phenomena, such as color-screening due to quark pair 
creation and many-body confinement forces are expected to play an important role. These 
long-range effects are difficult to accomodate within the standard two-body confinement force 
model of Lipkin-type with a universal confinement strength for baryons and compact six- 
quark states. In order to model the complex color dynamics in a compound six-quark system, 
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we have proposed that the effective two-body confining strength in the d! is weaker than in a 
single nucleon (model III). We then obtain a d' mass, = 2092 MeV, which is compatible 
with the experimentally suggested resonance mass. In the absence of a solvable theory of 
confinement, our harmonic oscillator model III should be viewed as an attempt to model 
the more complicated color dynamics in a genuine six-quark system. Our conjecture of a 
weaker confinement strength is based on the observation that the average distance between 
any pair of quarks in the d 1 is significantly larger than in the nucleon. In this long-range 
regime the confinement interaction is poorly understood. A harmonic oscillator confinement 
model clearly shows the inverse relation between the size of the system and the confinement 
strength. If the d' is experimentally confirmed, it would be an indication that the effective 
two-body confinement strength is weaker in a genuine six-quark system with the size of the 
d'. 

It remains to be seen whether our hypothesis of a weaker confinement strength in a 
compound six-quark system affects other results in the B = 2 sector such as NN scattering 



phase shifts or the deuteron electromagnetic form factors [49 . According to our experience, 



these observables are rather insensitive to the details of confinement, and we do not expect 
any qualitative changes of previous results. Nevertheless, new calculations, e.g. of deuteron 
observables should be carried out to explicitly check the consequences of our conjecture. 

An important test of the d 1 hypothesis is the simultaneous description of the mass and 
the free pionic decay width ~ 0.5MeV of the d' with the same set of parameters. Model 
III gives a free d 1 decay width of IV « 0.3 MeV, which is in agreement with experimental 
result ||28|| . In the future, a detailed analysis of the Nd' and d'd' interactions similar to the 
calculation of the NH 53] and HH |57j] interactions, as well as a calculation of the cross 
section for d'N — > NNN should be done. This will provide a stringent test of whether the 
theoretical d' radius 7> = 1.53 fm is consistent with the empirical 'in medium' d' decay 
width T medium « 5 MeV. 

On the experimental side, the search for a clear d! signature in elementary reactions, 
such as pp — > pp7r + 7T~ [[Tj|] and 7<i — ► pun []5"8|] , has begun. The experimental search for 



dibaryons provides a unique chance for probing the largely unexplored and little understood 
phenomenon of color confinement at larger interquark distances. In this regime our tradi- 
tional confinement models are likely to break down. If the existence of narrow dibaryons 
is confirmed, it will also have important implications for deep-inelastic electron scattering 



off nuclei [59], for the equation of state of nuclear matter at higher densities, and the inner 
structure of neutron stars [B^ . 
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Appendix 

In this appendix, we present all norm and Hamiltonian kernel matrix elements neccessary 
to diagonalize the algebraic generalized eigenvalue problem of Eq. (K% . We calculate the 
Hamiltonian and norm kernels 

HijisuSj) = (GCMilAHlGCMj) ; A/^s;, Sj) = (GCM i |.A|GCM i ), (31) 

where the abbreviation |GCMj) introduced in Eq. (|TjJ) denotes 

IGCM,) = -L[(lM) 3 / 4 ] 6 / dh ■ \z 8 ^> T =»> c = a x Y L=1 (s t )] J=0,T= ° ,C 

V47T J L 

x {nexp(-^(r.-|5)>)}{n«p(->, + |?)'; 

[ $ 5 T =1,T T =0 x [2n] C T =l $ 5 D =0,T C =0 x [ 11 ]J D =1] S = 1 ' T =° [222] C=0 ( g 2) 



X 



/S=1,T=0,C=0 



In Eq. (|52]) 2'5=i,t=o,c*=o j g S pi n (S)-isospin(T)-color(C) wave function of the d', where the 
diquark (Sd = 0, T D = 0, Cp = 1) and the tetraquark (St — 1, T T — Ct = 1) STC-wave 
functions are coupled to the total STC wave function of the d', with S = 1, T = 0, and total 
color C = of the color singlet state [222] i. The projection onto good angular momentum 
L = 1 of the orbital part of the matrix elements of Eq. ( PTj ) is done after the integration over 
all quark coordinates shown in Fig. £|: rit=i / dvi. Introducing the notation A = l/(26§), we 
obtain for the orbital integrals before projection, first for the norm matrix kernels Mif. 



■A/y = ex P ~ 0* ~ s i' 2 



2A 
3~ 

N$ = exp (-^-( s i - s ?)) ex P (+y s * • s j) 



JS/J* = exp (- — (s 2 - s. 2 ) exp -— s, • Sj (33) 



— (s 2 - s 2 ) j exp (^- T 

Here and in the following, D refers to the direct kernel, X represents the one-quark-exchange 
diagrams for quark 4 and 6, and XX refers to the two-quark-exchange diagrams of the 
pairs (34)^>(56) in orbital space. The subscripts in the kinetic energy and the two-body 
interaction kernels denote the quark coordinate. 

2 

The non-projected kernels of the kinetic energy operator Tj = = — ^~V 2 read 
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(34) 



where the norm kernels Mij of Eq. fl33|) have been factored out. For the two-body potential 
kernels depicted in Fig. [5] (for the case of the one-gluon-exchange potential), we arrive after 
some straightforward algebra at: 

/ A \ 3/2 ^oo 
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(35) 
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x £(-)' • / rf " 2 exp(-Ar 2 ) \/(r)^(^ Sl )^(^ Sj ; 

7 ™ " 



/.//; 

v£(*n* i ) = {N$IN§)-vj>{* u * j ) 



/ /I \ 3/2 ^oo 

\/ 1 ^(s,,s J )=^ X («^,) = ^6 X (^^)=^(-j (4tt) jf dr r 2 exp(-Ar 2 ) V{r) 

,j) = M xx (^J /2 (4vr) exp (- jsfj J™ dr r 2 exp(-Ar 2 ) V(r)i (Ar Sj ) 

A 



^23 ( S *5 S 



/ A \ 3 / 2 / A \ r°° 

V xx {s h s j )=N? x (^) (4tt) exp (--s 2 J jf drr 2 exp(-Ar 2 ) V{r)i {Ar Sl 
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(36) 
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In order to project all kernels onto good angular momentum L — 1 in the relative 
coordinate, we first expand the scalar products of the generator coordinates (sj-Sj) appearing 
in the exponential functions in terms of spherical harmonics 



exp {-aAsi ■ Sj) = 4tt ]T(-) ? i ii(aA Si Sj)Y l (si) ■ Y l (sj) 



(37) 



The projection may then be performed by application of the integral operator Pq C Ii = 

In the present case (L — 1, S — 1, T — 0) we obtain the following expressions for the 
projected norm kernels 
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(38) 



where the short-hand notation (STC) has been introduced to denote the reduced norm 
operator matrix elements (||1||), (| I-PJP"! |), and(| \Pi 



STC tdSTCi 



in spin-isospin-color space. 
Explicit expressions are give below. 

Some of the projected kinetic energy kernels are simply obtained by the substitution 
Afij — > Nfj OJ . The components proportional to ~ (sj • Sj) -Af^ are obtained after application 
of the trick: 
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(39) 

The replacement A/jj — > A/^ OJ applies also for most of the projected two-body potential 
matrix elements, except for 
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(40) 
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Finally, we explicitly give the spin-isospin-color (STC) factors for the different matrix 
elements, first for the norm kernel 

((S = 1,L = 1)J = 0,T = 1| 1 ||(S = 1,L = 1)J = 0,T = 0) = V3 
((S = 1,L = 1)J = 0,T = 0|Pf 6 TC ||(S = 1,L = 1)J = 0,T = 0) = 

((S = 1,L = 1)J = 0,T = 0|| Pg TC Pg TC ||(S = 1,L = 1)J = 0,T = 0) 

The same (STC) matrix elements apply to the kinetic energy and for the one-sigma-exchange 
potential, because as the norm operator, these operators do not depend on isospin, spin, or 
color. Table |V| collects all neccessary (STC) matrix elements for the remaining interactions. 
In the calculation of these matrix elements, we have made use of the well-known identities 

\ a i -\ a j = - 2 -+2P° , * i -* j = -l + 2P? j , T i -r j = -l + 2P? (42) 

and also of the fractional parentage decomposition of the tetraquark and diquark wave 
functions to evaluate the quark exchange operator acting in the different spaces. 
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TABLES 



Conf. 


Set 


b 3 [fm] 




ai 3) 







I 


0.603 


1.540 


24.94 MeV/fm 2 


= a< 3) 


Q 


II 


0.595 


0.958 


13.66 MeV/fm 2 




Q 


III 


0.595 


0.958 


13.66 MeV/fm 2 


5.00 MeV/fm 2 


L 


IV 


0.609 


1.060 


22.03 MeV/fm 


= <£> 


R 


V 


0.617 


1.122 


26.12 MeV/fm 2 / 3 


= af 


E 


VI 


0.648 


1.374 


46.70 MeV 





TABLE I. Quark model parameters. Sets I - III are obtained for the quadratic confinement 
potential (Q) of Eq. (7). Set I: without it- and cr-exchange potentials. Set II: with regularized ir- 
and c-meson exchange potentials. Set III: same as set II but with reduced confinement strength 
in the six-quark sector. Set IV: with linear confinement (L) of Eq. (8) instead of quadratic con- 
finement. Set V: with r 2//3 confinement interaction (R) of Eq. (9). Set VI: with the color-screened 

error-function confinement (E) of Eq. (10) using an inverse screening length of l//x= 0.8 fm. In 

(3) 

columns 5 and 6 we list the confinement strength a c as determined from the nucleon mass and 
the confinement strength employed for the calculation of the dibaryon mass and wave function 
as discussed in Section IV. 



Conf. 


Set 


M 2q 


M 4 g 


M d , 


b 6 








[MeV] 


[MeV] 


[MeV] 


[fm] 


[MeV] 


Q 


I 


636 


1501 


2610 


0.70 


2490 


Q 


II 


643 


1456 


2440 


0.75 


2316 


Q 


III 


621 


1309 


2092 


0.95 


2013 


L 


IV 


650 


1431 


2354 


0.86 


2261 


R 


V 


653 


1419 


2313 


0.93 


2234 


E 


VI 


661 


1422 


2288 


1.08 


2235 



TABLE II. The mass and characteristic size of the d! as obtained according to Eq.(18) 
for the six parameter sets of Table I. The diquark and tetraquark masses are also shown. The masses 
in the last column are obtained in a calculation without quark exchange between the diquark and 
the tetraquark clusters. As in table I, the abbreviations Q, L, R, and E identify the four different 
confinement models used in this work. The experimental mass of the d! is = 2065 ± 5 MeV [3]. 
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Conf. 


Set 




Vconf 


^gluon 






M d , 






[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 





I 


1048 


725 


-1041 






2610 






(997) 


(678) 


(-1063) 






(2490) 





II 


906 


460 


-617 


-44 


-143 


2440 






(866) 


(429) 


(-634) 


(-86) 


(-137) 


(2316) 


o 


III 


566 


270 


-517 


-23 


-82 


2092 






(545) 


(249) 


(-537) 


(-42) 


(-79) 


(2013) 


L 


IV 


677 


547 


-615 


-32 


-101 


2354 






(651) 


(524) 


(-636) 


(-57) 


(-98) 


(2261) 


R 


V 


574 


580 


-610 


-27 


-82 


2313 






(554) 


(562) 


(-634) 


(-45) 


(-81) 


(2234) 


E 


VI 


396 


722 


-639 


-21 


-49 


2288 






(402) 


(719) 


(-682) 


(-28) 


(-54) 


(2235) 



TABLE III. The kinetic (without rest masses) and potential energy contributions to the mass 
of the d' for the six parameter sets of Table I. The numbers in parentheses are obtained in a 



calculation without quark exchange between the diquark and the tetraquark clusters. 
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Conf. 


Set 


r D 




R RGM 






r RGM 
qq 


r TISM 
qq 






ffml 
[imj 


ffml 
[imj 


ffml 
[imj 


ffml 
[imj 


ffml 
[imj 


ffml 


ffml 
[imj 


Q 


I 


0.61 


0.74 


1.01 
(1.11) 


1.16 


0.96 


1.32 


1.38 


Q 


II 


0.65 


0.80 


1.10 
(1.21) 


1.24 


1.03 


1.42 


1.47 


Q 


III 


0.82 


1.01 


1.39 
(1.50) 


1.53 


1.30 


1.80 


1.82 


L 


IV 


0.74 


0.91 


1.32 
(1.43) 


1.41 


1.18 


1.67 


1.73 


R 


V 


0.81 


0.99 


1.46 
(1.56) 


1.53 


1.27 


1.84 


1.86 


E 


VI 


0.93 


1.15 


2.08 
(1.94) 


1.85 


1.48 


2.52 


2.25 



TABLE IV. The rms radii of the diquark (ro), tetraquark (ry), and the intercluster wave 
function (R§ GM ) with and without (in parantheses) the quark exchange diagrams. The total d! 
matter radius (see Eq.(21)) is also given. In the RGM formalism the same value for the oscillator 
parameter is used for the diquark, tetraquark, and relative motion wave functions. The mean 
distance between the colored clusters, R^ GM , is compared to the corresponding shell model radius 
RHO = ( 15 / 8 ) The mean distance between two quarks in the RGM (r^ M ) and TISM (rJ q ISM ) 
approach are also shown for comparison. The same notation as in Table I is used. 
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Diagram 


Operator Oij 


\ ■ Aj 


A? • A|<7i • a,- 




direct 


Cl2 


-4/V3 
-8/V3 
-2/V3 


8/(3^3) 
+8^3 



9^3 



lq-exchange 


n f>STC 

n f>STC 
^56 -^46 

^45 *4 6 
046 Pf 6 TC 
035 Pf 6 TC 


1/(2V3) 
-1/(2^3) 


-1/(2V3) 


V3/2 



-2/(3V3) 
2/(3 V3) 


2/(3 V3) 


3^/3/2 
4/(3^3) 


-5/(4V3) 
5/(4^3) 


5/(4>/3) 



5/(4^3) 


2q-exchange 


n f>STCf>STC 

V12 *46 F 35 

034 pf 6 Tc pf 5 TC 

05 6 Pf 6 TC Pf5 TC 
n pSTCpSTC* 

o, 5 P!l c P!F 


1/V3 
-2/V3 
-2/V3 
-5/(4V3) 
-5/(4>/3) 
1/(4V3) 
1/(4^3) 


1/V3 

2^3 

2^3 










-3^3/4 
9\/3/4 
9\/3/4 





-V3/4 
-V3/4 



TABLE V. Spin(S)-Isospin(T)-Color(C) matrix elements of the interaction kernels 
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FIGURES 



FIG. 1. The formation of a compound six-quark system from the initial ir + nn system. The d 
quark in the 7r + annihilates with a d quark in one of the neutrons leading to a q 2 ~q 4 clusterized 
six-quark state. In the stringlike bag model, colored diquark and tetraquark clusters are bound 
together by a color-electric flux-tube (string), which is rotating with angular momentum L = 1. In 
the present model, the relative q 2 — q A wave function is calculated from the Hamiltonian of Eq.(l). 



Dq 




D 



FIG. 2. The residual (a) one-gluon, (b) one-pion, and (c) one-sigma exchange interactions. The 
finite size of the quark-meson vertex described by the cut-off A is indicated by a small black dot. 
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FIG. 3. The radial dependence of typical confinement models in the nucleon. 




r [fm] 
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FIG. 5. The direct, one-quark Vij P2q STC and two-quark Vij P^ STC P§^ STC exchange diagrams 
for the one-gluon exchange potential. Corresponding diagrams are calculated for the confinement, 
one-pion, and one-sigma exchange potentials. The stringlike bag model includes only the diagrams 
in the first row. The quark exchange diagrams neccessary to satisfy the Pauli principle for the whole 
six-quark system are neglected in the stringlike bag model but included in the present theory. The 
quark exchange kernels increase the d! mass by about 100 MeV. 
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FIG. 6. The relative wave RGM function between the tetraquark and diquark clusters with 
(Pauli-on) and without (Pauli-off ) inclusion of the quark exchange diagrams for the smallest (model 
II) and largest (model VI) ci'-dibaryon. The RGM wave functions are compared to the pure 
harmonic oscillator wave function of Eq. (23). The upper graph shows the results for the quadratic 
confinement (model II) , while the lower graph is obtained for the color-screened confinement (model 
VI). For model II (upper graph) the clusters overlap, while they are just touching in model VI. 
The corresponding rms radii are given in Table IV. The wave functions of the other models are in 
between these two extremes. 
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